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I.  INTRODUCTION 


The  free  electron  laser  (FEL)  is  an  adaptable  source  of  coherent  radiation 
that  can  operate  at  high  power  and  has  many  potential  applications  in  scientific 
research,  industry,  national  defense  and  medicine.  The  FEL  uses  a  high  quality 
relativistic  electron  beam  made  to  oscillate  by  passing  through  a  periodic 
magnetic  field  to  amplify  coherent  optical  radiation  [Ref.  1].  The  basic  FEL 
configuration  (Figure  1)  consists  of  an  electron  accelerator  that  injects  a  high- 
quality  electron  beam  into  the  transverse  periodic  magnetic  field  called  an 
"undulator."  In  the  undulator,  the  interaction  between  the  beam  and  the  co¬ 
propagating  radiation  takes  place.  The  electrons  are  the  power  source  in  the 
FEL,  and  can  provide  more  than  a  gigawatt  of  peak  power  with  an  average 
power  ranging  from  kilowatts  to  ten  megawatts.  The  FEL  converts  the  electron 
beam  kinetic  energy  to  electromagnetic  radiation,  and  can  be  made  to  operate 
as  an  amplifier,  or  as  an  oscillator.  In  the  first  case,  a  high-current  beam 
amplifies  a  low-power  optical  wave  during  a  single  pass  through  the  interaction 
region  without  a  need  for  a  optical  feedback  system.  In  the  second  case,  mirrors 
placed  beyond  each  end  of  the  undulator  form  a  resonator,  or  optical  cavity.  In 
this  cavity,  radiation  is  reflected  between  the  mirrors.  Oscillating  electrons  are 
injected  into  the  cavity  to  overlap  the  rebounding  optical  pulse.  The  laser  field 
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grows  on  each  pass,  and  becomes  eventually  large.  In  both  cases,  H  is  critical 
that  the  electrons  spatially  overlap  the  optical  mode  for  energy  exchange  to 
occur  between  the  electron  beam  and  the  electric  field  of  the  light  wave  [Refs. 
2.3]. 


beomdump,  or  reoccelerotor 


radiation 

leaves 


Figure  1.  Schematic  of  the  Interacting  Elements  of  the  FEL  Oscillator. 
Successive  electron  pulses  enter  and  propagate  through  the  transverse 
periodic  magnetic  field  set  up  by  the  unduiator.  The  superimposed 
rebounding  coherent  optical  wave  is  amplified  by  the  electron  beam  energy 
exchange. 


The  basic  physics  of  the  FEL  is  analyzed  theoretically  by  solving  the 
relativistic  electron  particle  dynamics,  and  coupling  the  solutions  to  Maxwell’s 
optical  wave  equation  [Ref.  4,5].  The  three  interacting  elements  of  a  high-quality 
relativistic  electron  beam,  undulator-magnetic  field  that  causes  electrons  to 
wiggle,  and  resonant  optical  cavity  to  provide  feedback,  combine  to  produce  light 
at  wavelengths  that  range  from  the  microwave  to  the  ultraviolet.  The  forward 
propagating  electromagnetic  wave  and  a  transverse  magnetic  field,  whose 
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electric  and  magnetic  fields  are  perpendicular  to  the  direction  of  propagation, 
give  rise  to  an  axial  force  that  can  extract  energy  from  the  electron  beam. 

The  basic  mechanism  of  the  coherent  energy  exchange  is  the  bunching  of 
the  electrons  at  optical  wavelengths.  The  electrons  leaving  the  accelerator  are 
randomly  positioned  over  many  optical  wavelengths.  There  are  typically 
approximately  10^  electrons,  or  more,  in  each  section  of  the  electron  beam  one 
optical  wavelength  of  light  long.  As  the  light  and  electrons  interact  at  the 
beginning  of  the  undulator,  some  electrons  gain  energy  and  some  lose  energy. 
Those  that  gain  energy  move  a  little  faster  longitudinally  and  those  that  lose 
energy  move  a  little  slower;  this  eventually  creates  a  bunch  in  each  optical 
wavelength.  See  Figure  2.  The  bunching  process,  or  localized  spatial 
positioning  of  the  electrons  within  a  optical  wavelength,  is  responsible  for 
coherent  radiation  at  the  end  of  the  undulator  and  creates  the  feedback  process 
[Ref.  5]. 

As  a  consequence  of  the  non-linear  properties  of  the  gain  mechanism, 
higher  harmonics  of  electromagnetic  waves  can  be  generated  in  FEL  oscillators 
and  amplifiers  [Ref.  6,7].  The  generation  of  harmonics  in  FELs  can  be  used  to 
reach  optical  wavlengths  with  moderate  to  low  electron  beam  energy.  This  work 
discusses  the  advantages  of  using  harmonics  and  contributes  to  the 
development  of  compact  FELs  where  size  and  complexity  can  be  substantially 
reduced.  The  compact  FEL  would  use  a  micro-undulator  with  a  length  in  the 
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I  Optical  Wavelength 


Beginning  of  Undulator 


End  of  Undulator 


Figure  2.  Bunching  Process  of  Electrons  in  the  PEL 
Some  of  the  electrons  gain  energy  and  some  lose  energy  within  one  optical 
wavelength  as  they  travel  through  the  undulator.  The  electrons  spatially 
overlap  due  to  their  different  energies  and  cause  amplification  of  light  at  the 
end  of  the  undulator. 


range  of  ten  centimeters.  The  compact  FEL  would  subsequently  be  made 
availabe  to  a  much  broader  range  of  users  [Ref.  8].  The  considerable  effort 
directed  to  the  development  of  this  device  is  motivated  by  the  need  to  reach 
short  wavelengths  in  the  most  inexpensive  manner. 

In  Chapter  II  the  development  of  the  mathematics  used  to  model  the  FEL 
interaction  and  describe  the  gain  mechanism  is  outlined  to  provide  the  necessary 
background  to  understand  the  basic  physics  of  the  FEL.  Chapter  III  discusses 
the  concept  of  the  compact  FEL  from  the  perspective  of  the  accelerator  and 
micro-undulator  components  of  the  FEL  system  and  motivates  the  realization  of 
these  devices.  The  evaluation  of  three  micro-undulators  operating  on  harmonics 
at  an  optical  wavelength  of  0.4  microns  are  next  explored  in  Chapter  IV.  An 
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analyses  of  the  gain  for  these  FEL  systems  and  a  careful  discussion  of  the 
dependence  on  the  electron  beam  quality  is  done  in  order  to  more  fully 
understand  the  physics  and  define  new  limits  on  various  parameters. 

The  author’s  most  significant  contributions  made  toward  the  operational 
advancement  of  a  Compact  Free  Electron  Laser  are  summarized  below. 


•  Micro-undulator  operation  for  proposed  compact  FEL  devices  are  simulated 
to  explore  the  wider  use  of  harmonics  in  order  to  reach  optical  wavelengths 
with  moderate-to-low  electron  beam  energy. 

•  Gain  degradation  due  to  energy  spread  and  strong  field  saturation  effects 
are  shown  to  demonstrate  the  advantages  of  using  harmonics  in  order  to 
reduce  size  and  complexity  of  FEL  systems. 

•  The  Compact  FEL  designs  are  explored  by  five  methods  that  model  FEL 
physics  with  different  levels  of  sophistication  and  mathematical  formalism. 
The  merits  of  predicating  gain  by  each  method  are  identified  to  establish 
their  effectiveness  to  explore  new  designs. 

•  The  data  points  on  the  curves  in  Figures  1 5  through  20  represent  the  result 
of  several  computer  simulations  that  calculate  the  FEL  gain  using  either  a 
CRAY  computer  Sun  Workstation,  or  Personal  Computer. 

•  An  analysis  of  beam  quality  is  used  to  arrive  at  new  electron  beam  size 
limits  that  simplify  gain  analysis  for  all  undulator  designs. 
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II.  BASIC  PHYSICS  OF  THE  FREE  ELECTRON  LASER 


A.  DYNAMICS  OF  ELECTRON  TRAJECTORIES 

This  section  contains  the  derivation  of  the  electron  pendulum  equation  [Ref. 
2-6],  that  calculates  the  dynamics  of  individual  relativisic  electrons  as  they  are 
affected  by  the  electric  and  magnetic  fields  in  the  laser  cavity.  The  assumptions, 
approximations,  and  mathematical  steps  are  presented  in  detail.  Each  electron’s 
initial  conditions  determine  the  evolution  of  its  velocity  and  position.  This  motion 
can  be  represented  in  periodic  pendulum  phase-space  by  a  section  of  phase- 
space  one  optical  wavelength  long.  The  particular  path  each  electron  takes  is 
determined  by  the  amplitude  of  the  undulator  field  and  the  optical  fields.  The 
physics  can  be  understood  by  appealing  to  electron  phase-space  diagrams  that 
consider  the  state  of  electron  motion  in  the  two-dimensional  space  of  position  and 
velocity.  These  diagrams  will  be  utilized  to  study  the  FEL  interaction  and  provide 
a  method  of  evaluating  the  gain  mechanism  in  the  compact  designs. 

The  electron  evolution  through  an  operating  FEL’s  periodic  undulator  in  cgs 
units,  are  governed  by  the  Lorentz  force  equation: 

4(yP)  =  — W 

at  me 

where  e  =  magnitude  of  electron  charge, 
m  =  mass  of  electron, 
c  =  speed  of  light. 
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gj,  =  electron  velocity. 

B  =  optical  magnetic  field, 

E  =  optical  electric  field, 

=  external  magnetic  field  in  undulator, 

Y  =  dimensionless  Lorentz  factor  =  p  _  ^  ^  j-1/2^ 

The  undulator  field  seen  by  electrons  can  have  circular  polarization  where  the 
electron  path  in  the  laboratory  frame  of  reference  is  helical  along  the  z-axis.  An 
alternative  is  an  undulator  with  linear  polarization  where  the  electron  path  is 
sinusoidal  in  the  x-z  plane.  See  Figure  3. 


Electron  Pofh  in  Helical  Electron  Poth  in  Lineor 

Undulotor  Undulotor 


Figure  3.  Schematics  for  Perfect  Electron  Trajectories  in  Circular  and 

Linear  Polarized  Unduiators. 


The  ideal  helical  undulator  field  near  the  z*axis  in  rectangular  coordinates  is 

fi.  =  I l  =  fi  I COSM,  sinM,  0 )  (2) 

where  B  =  undulator  peak  magnetic  field  amplitude, 
kg  ==  undulator  wave  number  =  2-kIX^, 

=  undulator  period, 

z  =  longitudinal  position  of  electron  along  undulator. 


It  is  important  to  realize  that  this  representation  of  is  accurate  only  near 
the  z-axis,  since  away  from  the  axis  the  transverse  must  bend  to  satisfy 

Maxwell’s  equation,  V  x  =  O. 

The  electric  and  magnetic  fields  of  the  radiation  in  the  cavity  are  assume  to 
be  a  plane,  circularly  polarized,  monochromatic  wave.  Both  these  fields  travel 
along  the  z-axis  with  the  electron.  The  fields  are  expressed  in  terms  of  the  optical 
peak  electric  field  amplitude  and  a  sinusoidal  disturbance: 

=£(r)  [cosijr, 0] 

=  £(t)l  Sin\{r,  COSij;,  0  J 


where 


'If 

k 

(i) 

<p(tj 


optical  peak  electric  field  amplitude, 
kz  -  oit  +  q)(t), 
optical  wave  number=2Tc/A.. 
carrier  frequency  Q=kc, 
optical  phase. 


The  simplest  and  most  fundamental  electromagnetic  waves  are  transverse,  plane 


waves  where  E  and  (p  are  independent  of  time  and  space  z.  The  goal  of  the 


calculation  is  to  determine  how  the  details  of  the  electron  trajectories  in  the 


undulator  are  coupled  to  the  light  wave  in  the  resonator  cavity  to  produce  a  large 
optical  power  output. 
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The  energy  change  for  the  relativistic  electron  in  a  radiation  field  is 

(4) 

at  me 

where  d(ymc^) I dt  is  the  rate  at  which  the  electron  loses  or  gains  energy.  An 
initial  distribution  of  the  electrons  and  velocities  at  the  entrance  of  the  undulator 
maybe  chosen  so  that  on  the  average  the  electron  beam  loses  energy  to  the 
radiation  field  in  one  pass  through  a  finite  length  of  undulator.  Figure  2,  in  the 
previous  chapter,  shows  a  distribution  of  electrons  on  the  optical  wavelength  scale 
where  an  alteration  of  z-velocities  provides  coherent  bunching  at  the  end  of  the 
undulator  to  provide  amplification  of  the  radiation  within  the  resonator. 

The  components  of  the  Lorentz  force  equation  and  electron  energy  equation 
are 


transverse  component  x-y  plane 

■7( =  —  [£(1  -  Pj)(cosnr,  -sin»|r.O)  *  p^( -sin(k^).cos(*^).0)] 
dt  me'-  ^ 


(5a) 


longitudinal  component  z- direction 

-X  ( yPJ  =  —  [£(P,cos4r  -  p  simjf)  ♦flCP^sinf*,^)  -  p  cos(i^)] 
at  me 


(5b) 


Y  =  —  £  I P,  cosijf  -  p  sirn|>  ] .  (®) 

me 

For  realitivistic  electrons,  Pj=1 ,  the  electric  and  magnetic  optical  fields  in  (5a) 
nearly  cancel.  This  implies  >  E(  1  -  p,)  when  P,  *  1 .  Therefore,  the  first  term 
of  the  transverse  component  of  the  Lorentz  force  may  be  neglected  for  relativistic 
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electrons.  The  transverse  electron  velocity  is  determined  almost  entirely  by  the 
static  magnetic  field.  This  assumption  allows  the  transverse  component  to  be 
approximated  by 


-7(yPJ  =  —  Bi  -sm(k^z),cos(k^z),0).  (7) 

at  me 

Exact  integration  of  this  equation  is  allowed  if  we  assume  perfect  injection  of  the 
electrons.  This  assumption  removes  the  constants  of  intergation  and  leads  to  the 
transverse  velocity; 


=  —  ( cos(  sin(  k^z),0)  (8) 

Y 

where  K  is  the  root-mean-square  (rms)  undulator  parameter  defined  by 


eBX^ 
K  = - ^ 


(9) 


where  B  =  B^f  ^f2  is  the  root-mean-square  undulator  magnetic  field  along  the  z- 
axis,  and  A,  is  the  undulator  period.  Perfectly  injected  electrons  enter  the 
undulator  exactly  on-axis  without  any  random  spread  in  angles  or  positions,  i.e., 
(x,y)  =  (0,0) ,  where  x  is  the  wiggle  direction  and  /  the  direction  of  the  magnetic  field. 

The  undulator  field  gives  the  electron  a  small  but  important  transverse 
"wiggle"  velocity  p^.  The  radiation  fields  alone  have  no  significant  effect  on  the 
electron’s  transverse  trajectory,  since  forces  due  to  Ej.  and  nearly  cancel  for 
relativistic  electrons,  y>‘'-  Since  total  energy  of  the  wave-particle  must  be 
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conserved,  energy  lost  by  the  electron  is  gained  by  the  light  wave  passing  through 
the  cavity  according  to  (4).  The  relative  orientation  between  p  and  is 
determined  by  the  electron’s  initial  position,  within  a  radiation  wavelength,  A. 
The  determination  of  whether  the  electron  gains  or  loses  energy  depends  on  the 
sign  of  Y  in  (4).  When  y  is  positive,  the  electron  gains  energy  and  radiation  is 
absorbed.  When  y  is  negative,  the  electron  loses  energy  as  the  radiation  field 
grows,  causing  stimulated  emission.  The  electron  transverse  velocity  (8),  and  the 
electric  field  =  cosi);,  -sinijr ,  0 )  with  a z-component  can  be  used  to  solve  the 
energy  equation  for  y(t).  Substituting  into  (4),  the  expressions  for  andp^ 
yields 

Y  =  — [£(cos4r,  -sintif.O)]  •[(  — ]cos(*^z),sin(*^z),0  . 
me  11  Y  i 

Performing  the  vector  dot  product  operation,  and  making  use  of  the  trigonometric 
identity  cos(  ^  -  cos\|;  cos( )  -  sinij;  sin( k^z )  the  above  equation  can  be 

written  as 

Y  =  OOS(  c  +  <P  ) 
yme 

where  li=(k  kj  z  -  o)t  \s  the  electron  phase,  and  k^  is  the  phase  of  the  optical 
wave.  Differentiating  C  with  respect  to  z,  for  any  f>0.  we  see  that 
AC  ={f<^k^)  Az.  For  relativistic  electrons /3^=  7,  the  optical  wavelength  A  is  much 
smaller  than  the  undulator  period  A,,  implying  k»k^.  Therefore,  AC  ~  k^z 
describes  the  position  of  electrons  on  the  optical  wavelength  scale  with  respect  to 
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the  combined  optical  and  undulator  forces.  The  actual  electron  beam  initial 
conditions  can  be  chosen  so  that  the  beam  predominantly  loses  energy  to  the 
radiation  field  as  the  electrons  travel  through  a  finite  length  undulator.  This  gain 
mechanism  will  be  discussed  in  the  next  section. 

A  more  useful  form  of  (10)  can  be  obtained  by  relating  the  electron’s  energy 
to  the  electron’s  position  in  a  section  of  the  electron’s  beam  one  optical 
wavelength  long.  Equation  (8)  and  the  definition  of  the  Lorentz  factor  in  (1)  can 
be  equivalently  written  as  y  =  1  ~  Pf  "  p!  where  !  y^.  By  eliminating 

P!  =  and  solving  for  p/,  we  have 


P?  =  1  -  - 


(11) 


Taking  the  differential  of  with  respect  to  time,  y  is  related  to  p^  by  the  following 
relationship 


Y  1 


P.P.- 


(12) 


In  order  to  relate  y  to  C .  we  must  take  the  differential  of  C  =  ( it  +  )  z  -  of . 
the  electron  phase  velocity  with  respect  to  time,  and  solve  for  p^.  The  rate  of 
change  of  the  electron’s  velocity  in  the  z-direction  is  given  by  p^.  Substituting  this 
result  into  (12).  we  see  that 


Y  ^  '  {k^kjc 

Recalling  k»k^  and  /3^~7.  (13)  can  be  written  as 


(13) 
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(14) 


Y  ■  C 

Y  (1  *K^)  w 

where  w=/fc. 

At  resonance,  the  electron  has  zero  phase  velocity  v=0  and  will  pass  through 
one  undulator  period  as  one  optical  wavelength  passes  over  it.  This  resonance 
condition  provides  the  maximum  coupling  and  energy  exchange  between  the 
electron's  transverse  velocity  and  the  electric  field  but  not  necessarily  the  maximum 
transfer  of  power  for  the  entire  beam.  See  Figure  4. 


Figure  4.  Resonance  Condition  in  the  Free  Electron  Laser. 

In  order  for  the  work  pj.  •  E,  to  be  positive  the  electron  must  move  back  an 
optical  wavelength  in  transversing  one  undulator  period. 


For  slow,  efficient  exchange  to  occur  between  the  electron  beam  and  the 
optical  wave  the  electron  phase  C  must  remain  close  to  0,  implying 
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( A:  +  )  p,c  -  (0  =0  or 


-pj. 


(15) 


By  making  use  of  binomial  expansion  for  y>1.  ('ll)  becomes 


K- 


1 


Y 


2=1  -I 


1 

\  Y^ 


or  (1  -PJ  =A 


1 

l^Y^ 


Substitution  of  this  expression  for  (1-p^)  into  (15)  yields 


1 

2y2 


(16) 


The  resonance  condition  in  terms  of  optical  and  undulator  wavelengths  for 
relativistic  electrons,  P^~7,  is 


A  =  A, 


'  1 


2y' 


(17) 


Finally  by  making  use  of  (10),  (14),  and  the  resonance  condition,  the  non¬ 
linear  pendulum  equation  describing  electron  particle  dynamics  in  the  FEL  takes 
on  the  form; 


C  = 


2€KE<jj„ 

yhnc 


cos(C 


<P) 


(18) 


where  =  k^c  =  2i:cl 

The  FEL  interaction  takes  place  while  the  light  and  electrons  travel  through 
the  undulator.  In  the  relativistic  limit,  the  interaction  time  is  close  to  Lie  for  both 
the  light  and  the  electron,  where  L  is  the  length  of  the  undulator.  Define  x=ct/L 
as  the  dimensionless  interaction  time  of  the  electron  with  the  optical  field  as  they 
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travel  through  the  length  of  the  undulator  where  0<  t  <  1 .  The  Pendulum  equation 
takes  the  form: 


COS  ( C  +  <p ) 


(19) 


where  H  is  the  number  of  undulator  periods. 
Let 


\a  = 


ATzNeKLE 
y^mc^  / 


.  C- 


dK. 


_  dC 


——and  V  =  C  = 
dx^  dx 


then  the  final  equation  becomes: 


Pendulum  Equation 


The  dimensionless  optical  field  amplitude  |a|  gives  the  strength  of  the 
coupling  that  occurs  between  C  and  v,  and  determines  the  electron  bunching  rate. 
For  each  electron  traveling  through  the  undulator,  resonance  is  measured  by  the 
electron  phase  velocity  v.  The  phase  velocity  is  the  difference  between  the 
frequency  at  which  electrons  pass  over  wavelengths  of  the  undulator  and 
frequency  at  which  wavelengths  of  light  pass  over  the  electrons  [Ref.  6],  The 
Pendulum  equation  describes  the  electron  evolution  where  the  initial  electron’s 
phase  and  phase  velocity  is  defined  by  C(t=0)  =  Co  and  v(t=0)  =  v^. 


B.  GAIN  MECHANISM  FOR  THE  LOW-CURRENT,  LOW-FIELD  FEL 

A  necessary  attribute  of  the  compact  FEL  is  the  ability  to  have  a  sufficient 
fractional  increase  in  the  optical  power  per  pass  through  the  undulator.  This 
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requirement  is  important  so  that  losses  due  to  the  low  reflectivity  and  other  optical 
resonator  losses  can  be  overcome  [Ref.  9].  The  laser  oscillator  is  designed  with 
mirrors  carefully  aligned  so  that  light  can  bounce  back  and  forth  between  these 
mirrors  and  have  very  little  loss  per  bounce.  A  single  pass  of  the  electron  beam 
is  defined  as  the  transport  of  the  beam  through  the  optical  cavity  from  beginning 
to  end  of  the  undulator.  If  the  net  amplification  between  the  mirrors,  including 
scattering  losses,  exceeds  the  net  reflection  loss  at  the  mirrors,  the  coherent 
optical  power  will  build  up  in  the  system.  A  measure  of  the  fractional  increase  in 
the  optical  power  is  called  the  gain  G  and  defined  by 


Opt 


where  =  change  of  power  along  the  undulator 

Pop,  =  power  at  beginning  of  undulator. 

The  gain  can  also  be  expressed  as  the  ratio  of  the  energy  lost  by  the  electron 

beam  to  the  power  required  to  create  the  energy  lost  by  the  electron  beam.  It  is 

assumed  that  the  electrons  do  not  interact  with  any  other  medium  so  that  when 

the  electron  beam  loses  energy,  all  the  energy  is  transferred  to  the  light  wave. 

This  is  critical  in  the  start-up  of  the  PEL  where  the  strength  of  the  optical  field  is 

initially  very  small.  The  following  discussion  highlights  the  derivation,  assumptions 

and  approximations  of  the  low  current-weak  field  gain  formula  [Ref.  4,5]. 

When  the  optical  field  is  considered  weak,  a«-K,  there  is  little  change  in  the 
electron  phase  velocity  during  the  PEL  interaction.  On  the  first  few  passes  of  light 
in  the  oscillator,  the  optical  field  strength  |e|,  starts  out  very  small.  When  the 
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dimensionless  current  is  small,  y«1,  there  is  little  change  in  the  optical  field 
amplitude  or  phase  in  a  single  pass. 

The  assumptions  used  in  this  calculation  of  gain  are: 

(1)  the  PEL  is  already  operating  and 

(2)  the  gain  is  low  so  that  a  -  small  - 0. 

o  o 

Using  the  pendulum  equation  C  =  v  =  |a|  cos  (  C  *>(>),  Ref-  5  solves  for  v(i),  the 
electron  phase  velocity,  by  using  perturbation  theory.  This  expression  is  given  by 


V  (t)  =  *  ^[sln(C,  +  -  sin(C„)]  - 

i 

^  -  i(cos{2C,  *2v,t)  -cos(2<:j)  *cos(v,t)  -1  -  v^TSin{ C J C50s( *  v„t ) 

\i'*  ^ 


where  the  initial  optical  field  is  \a(o)  \  =a^  and  (p(o)  =0. 

The  net  change  in  energy  of  the  electron  beam  is  defined  as  "  the  energy 
electrons  gain  as  work  is  done  on  them  by  the  radiation  field"  minus  "the  energy 
electrons  lose  to  the  radiation  field."  A  net  change  in  energy  of  the  electron  beam 
does  not  occur  for  the  term  proportional  to  because  just  as  many  electrons  lose 
energy  as  gain  energy.  The  higher-order  terms  proportional  to  in  (21)  will  shift 
and  skew  the  distribution  of  electron  phase  velocities  so  there  can  be  a  net 
change  in  the  energy  of  the  beam.  To  second  order  in  a^,  there  exists  a  net 
change  in  the  electron  energy  given  by 

2 

<v>  =  +  — —  [  2  cos  (  V^T  )  -  2  +  V^X  sin(  v^x  ;  ] .  (22) 

2vr 

where  <...>  is  used  to  denote  a  normalized  average  over  the  electrons. 
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The  average  energy  lost  by  an  electron  is  ymc^(<v>  -v^)IAnN  [Ref.  5] 
where  the  electrons  are  assumed  initially  distributed  uniformly  in  phase  and  mono- 
energetic  with  phase  velocity  v^.  The  energy  loss  by  all  electrons  in  the  beam 
of  volume  element  SV  is  therefore 


y  mc^  {<v>  -  v^) 
47:  iV 


(23) 


where  p.  is  the  electron  density.  For  low-gain  the  electrons  are  assumed  to 
remain  well  inside  the  optical  mode,  and  do  not  disturb  the  shape  of  the  optical 
beam.  The  electron  beam  energy  in  volume  element  SV  is  deposited  into  the 
optical  field  during  the  FEL  interaction.  The  radiation  energy  density  that  enters 
the  optical  field  from  the  electron  beam  is  governed  by 

E^bV  (24) 


5 


optical 


47t 


Where  =  radiation  energy  density 

E  =  magnitude  of  circularly  polarized  radiation  field 

bV  =  volume  element  of  optical  beam. 

The  gain  is  then  found  by  taking  the  ratio  of  (23)  and  (24),  and  equating  the 

energy  loss  by  the  electron  beam  to  the  energy  gained  by  the  light  wave.  The 

gain  is 


Substituting  (22)  into  the  above  expression,  and  defining  a  dimensionless  current 


~  2  2 
•J 


j  =6N {enKL)^pJ  y^mc^  the  gain  that  occurs  in  the  FEL  is 


(26) 


.  [2 -2cos(v„t) -v^TSin(v^T)| 

G=; - 5 - 

This  equation  holds  when  gain  is  low.  The  gain  spectrum  is  a  plot  of  the  final  gain 
at  T  =  1  verses  the  initial  phase  velocity  v^.  The  v^-axis  can  be  considered  as  a 
function  of  the  optical  wavelength  centered  at  the  resonant  wavelength  expressed 
by  (17).  The  gain  spectrum  is  anti-symmetric  in  for  j«^  with  a  peak  gain  of 

G  =  0.1357  (27) 

at  Vo=2.6.  At  resonance,  there  is  no  gain,  while  at  values  of  v<0  there  is  net 
absorption  of  the  optical  power.  A  change  in  resonance  of  Av=7i  can  shift  the 
interaction  from  amplification  to  absorption.  Shown  in  Figure  5  is  the  gain 
spectrum  for  moderate  current, /=5,  in  weak  optical  fields  3^=1.  The  gain  for  this 
design  is  80%.  The  gain  spectrum  is  no  longer  anti-symmetric  in  Vo  but  shifted  to 
the  left  where  some  gain  at  resonance  now  occurs.  Since  this  dimensionless 
current  corresponds  to  a  moderate  electron  beam  current,  deviations  of  the 
spectrum  from  the  low  current  case  are  expected. 

Two  corrections  to  (26)  include  (i)  a  coupling  factor  that  incorporates  Bessel 
functions  that  depend  on  K  and  the  undulator  polarization  and  (ii)  a  filling  factor 
that  includes  the  optical  mode  distortion  effects  [Ref.  5].  Both  these  factors  are 
significant  in  FEL  gain  calculations  and  will  be  explored  in  Chapter  IV. 

C.  DERIVATION  OF  THE  OPTICAL  WAVE  EQUATION 

We  can  derive  a  wave  equation  to  describe  the  FEL’s  optical  evolution.  The 
assumptions,  approximations  and  calculations  found  in  [Ref.  2-6]  are  reviewed. 
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****  Compact  FEL  Gain  Curve  **** 

j=5  a  =1  N=49 

o 

0.8 

0.0 

-0.8 

Figure  5.  Weak*Field  Gain  Spectrum  G(vJ  For  Moderate  Current 

The  purpose  of  the  review  is  to  provide  the  reader  a  basic  understanding  of  the 
evolution  of  light  in  the  presence  of  relativistic  electrons  that  travel  through  a 
magnetic  field  set  up  by  an  undulator. 

As  relativistic  electrons  initially  enter  the  undulator  field,  approximately  one 
photon  will  be  emitted  for  each  electron  that  passes  through  the  undulator.  This 
is  spontaneous  emission.  The  total  field  from  many  electrons  rapidly  establishes 
a  classical  light  wave.  The  number  of  electrons  contained  within  the  optical  volume 
element  5V  depends  on  the  electron  density,  and  the  ratio  of  the  electron  beam 
size  to  the  optical  mode  size.  As  further  amplification  of  the  light  wave  occurs, 
approximately  1 0®  photons  per  pass  are  emitted  into  the  optical  volume  element 
causing  the  light  to  become  a  coherent  emission  source.  This  electromagnetic 
radiation  provides  the  driving  force  on  the  oscillating  electrons  as  they  travel 
through  the  undulator.  Using  Maxwell’s  equations,  the  evolution  of  the  amplitude 
and  phase  of  the  coherent  optical  wave  can  be  coupled  to  the  single-particle 
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Lorentz  force  equations.  The  electric  field  E(z,t)  changes  in  position  and  time  due 
to  the  amplification  process  and  the  waveform  passing  over  the  electrons.  The 
electric  component  of  the  optical  field  (3)  can  be  approximately  expressed  by  a 
radiation  vector  potential  A  : 

A  =  ( sin^r,  COS^r,  0 ) ,  (28) 

k 

where  E(t)  =  electric  field  vector, 
k  =  2Tz/k, 

ilr  =  kz  -  o)t  +  <p(t)  and  <p(t)  =  optical  phase. 

The  amplitude  of  the  electric  field  is  assumed  to  be  independent  of  the  position  z. 
It  is  also  assumed  that  the  electrons  remain  well  inside  the  optical  mode  waist,  and 
therefore,  no  transverse  dimension  considerations  are  included  in  the  radiation 
vector  potential. 

By  making  use  of  the  Lorentz  gauge,  the  inhomogeneous  wave  equation 
[Ref.  10]  becomes 

(29) 

[  a®  J  c  * 

Where  -  {M c^)  ^1  dt^]  =  wave  operator, 

c^/dy\  e^/dz^) 

=  transverse  election  current  density  source  term. 

This  expression  describes  a  plane  wave  traveling  in  the  z-direction  when  the 
amplitude  and  phase  are  held  fix.  When  the  amplitude  and  phase  slowly  evolve 
over  an  optical  wavelength,  first  and  second  derivatives  in  the  spatial  and  time 
domains,  and  terms  containing  two  derivatives  can  be  negle'-ted,  since  the 
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radiation  bandwidth  is  emitted  in  a  narrow  spectrum.  Substituting  (28)  into  (29) 
and  using  this  approximation,  the  left  side  of  the  wave  equation  can  be  simplified 
to  [Ref.  5] 

—  (cosiv, -siniir.O) ^  (siniif.costir.O).  (30) 
[  ar^  J  c  dt  c  dt 

Equations  that  are  slowly-varying  can  be  constructed  by  projecting  the  wave 
equation  on  two  unit  vectors  in  the  following  manner; 


where 


1  8^  I*  / .  V  ^  4^ 


dr 


€2  = - JJt)  eg 

c 


(31) 

(32) 


=  ( coS(}f,  -slrui;,  0 )  (33a) 

€2  =  ( sin([f,  cosijr,  0 ) .  (33b) 

The  fast  factors  'bos  (4r)"and  'fein  (4r)"each  oscillate  once  as  an  electron  near 


resonance  passes  through  an  undulator  period. 

The  substitution  of  (30)  for  [  -  ( 1  / )  ^ldt^]A(t)  and  (33a)  for  e,  are 

made  into  the  left  side  of  (31).  A  similar  procedure  is  performed  with  (33b)  for  62 
and  the  left  side  of  (32).  Performing  the  vector  dot  product  operation  on  the  left 
side,  Maxwell’s  equations  become: 
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(34) 


c  dt 


(35) 

c  dt  c  "  ^ 

The  source  term  J  {t)  describes  the  electron  current  density  as  a  function 
of  the  time  the  electrons  interact  with  the  optical  wave  and  the  longitudinal  position 


of  the  electron  in  the  undulator.  The  single  particle  current  density  is  determined 
by  [Ref.  10]; 


(36) 


where 


=  transverse  electron  current  density  source  term  of  the 
ith  electron 

e  =  charge  of  electron 

c  =  speed  of  light 

=  transverse  electron  velocity 

6^^  (Jt -fj  =  delta-function  that  describes  the  electron’s  location  at 
position  in  3-dimensional  space 
=  trajectory  of  the  ith  electron. 


Since  Maxwell’s  equations  are  driven  by  the  total  beam  current,  the  sum  of  all  the 


electrons  particle  currents  must  be  determined.  Equation  (36)  is  modified 


accordingly: 


I  I 

where  the  sum  is  over  all  i  electrons  located  in  a  small  volume  element  bV. 

The  electron’s  transverse  motion  is  almost  entirely  determined  by  the 
undulator  field,  so  (8)  can  accurately  be  used  to  express  in  the  total  beam  current. 
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Using  (8),  (36),  (33a)  and  (33b).  the  right  side  of  Maxwell’s  first-order  scaler 
equation,  (34),  is  transformed  into 

—  yj  r )  •  e,  =  (f  -  f.)  [cos(  k^z )  cosiii  -  sin(  k^z )  simi;  ] . 

^  Y  1 

Using  the  trigonometric  identity  cos{k^z  -•-Tlr)  =cos(i^z)  cosiir  -  sin  (i^)  sin  4r.  the 
final  form  of  the  transverse  electron  current  density  projected  along  786,  is 

/(f)!,  =  — 536<3)(i--f;)COS(A:^  +  ilr).  (38) 

Y  1 

Following  similar  manipulations  with  (8).  the  final  *orm  of  the  transverse  electron 
current  density  projected  along  eg  is 

J  (,)  .|  =  ££« £  - f,) sin(*,2  -  4- ) .  (39) 

Y  . 

Substituting  the  electron  current  density  equations  (38)  and  (39)  into  Maxwell’s  two 
first-order  scaler  equations.  (34)  and  (35),  we  acquire: 

l^  =  22£*£6B'(r-r,)cos(*„:*4.)  («) 

c  dt  Y  . 

In  order  to  proceed  further,  the  volume  element  where  the  FEL  interaction  occurs 
must  be  defined.  This  differential  volume.  bV,  determines  how  the  appropriate 
approximations  of  the  electric  field’s  slowly-varying  time  and  space  coordinates  are 
made,  and  also  provides  a  means  to  sample  all  the  electrons  that  are  contained 
in  this  space.  The  electron  and  optical  beams  are  both  considered  cylindrical 
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where  natural  diffraction  of  the  light  is  ignored,  and  electrons  are  well  inside  the 
optical  mode  waist  See  Figure  6. 

Light  Bfoni 
Cltclrofi  Btom 

Optical  Field  Strueturf 

Figure  6.  Differential  Volume  Element  In  the  FEL 
Simplification  to  the  light  and  electron  beams  are  made  so  the  electron  beam 
remains  well  Inside  the  optical  beam  where  both  are  taken  to  be  cylindrical 
In  shape.  Also,  shown  are  several  optical  wavelengths  inside  the  optical  field 
contained  in  AV. 

The  volume  element  is  chosen  so  that  it  contains  several  wavelengths  of  light 
in  which  the  optical  field  structure  (E,  remains  essentially  constant.  An  averaging 
process  of  (40)  over  this  volume  element  for  the  optical  field  amplitude  and 
electron  particle  positions  must  be  performed  to  remove  6 -functions.  By  averaging 
over  an  element  6V  that  is  much  smaller  than  the  size  of  the  optical  field  envelope, 
all  the  field  sites  will  have  essentially  the  same  slopes  and  time  derivatives.  This 
results  in  the  left  side  of  (40)  remaining  unchanged.  The  right  side  also  must  be 
averaged  over  many  optical  wavelengths  where  all  the  6 -functions  that  exist  within 
the  volume  element  AV  must  be  summed.  Since  the  average  is  over  many  optical 
wavelengths,  a  large  number  of  electrons  are  either  bunched  or  randomly 
positioned.  The  electron  phase  C  will  provide  a  means  of  locating  electrons  within 
0  to  2ic  in  phase  space.  The  function  cos{k^  *  i|r(0) )  »  cos  [( it,  jt)z]  at  f=0  and 
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<p(z,fj=0  can  be  approximated  as  cos{kz^)  for  k»k^.  Recalling,  for  any  t>0  and 
k»k^,  »|r  follows  the  electron  position  in  a  section  of  the  beam  one  optical 
wavelength  long,  AC  =  klz.  Therefore,  the  position  of  the  electron  within  an  optical 
wavelength  is  simply  given  by  /^=kz. 

The  sum  ^  -  fj  cos(k^z  +  4r)  electrons  in  6V  can  be  replaced  by  an 

I 

average  over  sampled  electrons  within  the  optical  wavelength 
-  rj  <cos(C  +  4>)>  where  +  »}>  -  C +4>  and  <cos(C*4>)>  is  the 

i 

averaged  sum.  The  sum  of  the  6 -functions  for  all  electrons  in  the  volume  element 
6V  can  be  converted  to  an  electron  particle  density,  p,,  that  will  properly  weight  the 
average  sum  of  the  sampled  electrons  <cos  (C  +  <p)>  and  correctly  determine 
how  much  the  optical  wave  is  driven  by  all  electrons  in  6V.  The  form  of  the  field 
envelope  is  assumed  to  depend  on  both  z  and  t.  Combining  these  assumptions 
and  approximations  (40)  can  be  written  as 

^  =  l5££*p<cos(  ?-♦)>•  («») 

dt  Y 

The  evolution  of  the  optical  phase  is  determined  by  performing  a  similar  weighted 
average  over  sample  electrons  in  the  microscopic  volume  element: 

5^  -f;)sin(/:^  -  «|r)  -P,<sin(c  +  <p)> 

l 

where  now  (41)  is  written  as 
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(42b) 


-2kC€K  r 

£_^  = - p^<sin(  C  <P  )>  • 

at  Y 

These  two  non-linear  equations  are  coupled  through  the  electron  pendulum 
equation,  and  describe  the  long-term  evolution  of  the  amplitude  and  phase  of  the 


laser  field.  Using  the  dimensionless  parameters 
AnNeKLE  ct 


a  = 


2  2 


j.SN{enKL)‘-^ 


and 


y  me 


defining  a  =  |  a  |  e**,  (42a)  is  put  in  complex  form: 


da 

dx 


=  -j<e''^> . 


Optical  Wave  Equation 


The  utility  of  this  equation  in  expressing  gain  can  be  immediately  seen  by 
considering  values  of  C  that  either  maximize  or  minimize  the  real  part  of 
Examining  this  equation,  we  can  see  that  bunching  of  the  beam  at  C  =  rt  will  drive 
the  amplitude  of  the  optical  wave,  and  lead  to  gain.  Growth  of  the  wave  also 
occurs  with  increased  dimensionless  current  /  and  depends  on  the  electron 
distribution  <...>.  Electrons  bunched  at  n  provide  maximum  gain  where  as 
electrons  bunched  at  n/2  provide  minimum  gain. 

The  dimensionless  current /  provides  the  coupling  between  the  electron  beam 
and  the  light  wave.  The  response  of  the  optical  wave  to  bunching  in  the  beam  is 
also  determined  from  /.  By  comparing  the  magnitudes  of  the  dimensionless 
quantities  to  certain  threshold  values,  a  simple  identification  of  PEL  gain  regimes 
can  be  made  [Ref.  5].  When  y<1 ,  the  PEL  gain  is  low  with  non-exponential  growth. 
When  />1,  the  PEL  gain  is  high  with  exponential  growth.  Por  values  of  j  not  within 
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either  of  these  regimes,  the  gain  does  not  have  either  of  these  distictive  features. 
The  proposed  compact  PEL  design  has  dimesionless  current  j= 5.  Therefore,  the 
correct  theoretical  approach  to  analyzing  gain  for  this  device  must  be  carefully 
chosen.  A  discussion  of  the  methods  used  to  evaluate  gain  is  addressed  in 
Chapter  IV. 

D.  STRONG  OPTICAL  FIELDS  AND  SATURATION  EFFECTS 

The  Maxwell-Lorentz  theory  of  the  PEL  combines  Maxwell’s  wave  equation 
for  the  optical  wave  and  the  relativistic  Lorentz  force  equation  for  the  electrons  to 
describe  laser  dynamics  in  both  strong  and  weak  optical  fields.  The  extension  of 
this  theory  into  the  strong  field  regime  allows  the  use  of  the  optical  wave  and 
electron  pendulum  equations  to  accurately  study  the  electron  phase-space 
evolution  and  gain  in  the  PEL  [Ref.  5].  This  extension  is  allowed  because  both 
equations  assume  there  are  no  significant  changes  in  the  electron  energy  ymc^. 
This  implies  the  electron  energy  is  not  followed  self-consistently.  Strong  optical 
fields  can  be  reached  after  many  passes  in  a  low-gain  oscillator,  or  at  the  end  of 
a  high-gain  amplifier.  When  the  optical  field  in  either  the  low  or  high-gain  system 
becomes  sufficiently  strong,  the  initially  uniform  electron  beam  becomes  over¬ 
bunched  and  complicates  the  electron  phase  distribution  /( C.  v  ).  Overbunching 
is  the  process  where  electrons  become  overpopulated  at  the  end  of  the  laser  and 
do  not  efficiently  drive  the  optical  wave.  This  over-bunching  can  destroy  the 
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coherent  bunching  responsible  for  the  initial  gain,  and  indicate  the  onset  of 
saturation  and  the  reduction  of  gain. 

Strong  fields  begin  when  [a\=Tz.  Recalling  the  electron  pendulum  equation, 
V  =  C  =  |al  COS ( cos ( C  ■"  <!>)).  We  see  that  a  bunch  formed  at  the  proper  phase 
for  gain,  C  ♦  «p  =  n,  would  soon  evolve  by  AC  ~  -  i  a  |  =  -  ir  to  the  phase  for 
absorption  C  <p  =0.  The  electron  phase  velocity  is  proportional  to  the  rate  of 
energy  transferred  from  the  electron  bunch  and  determines  the  amount  of  power 
transferred  into  the  optical  wave.  Similarly  when  C +<!)=-«,  the  electron  phase 
velocity  v  and  rate  of  energy  transferred  from  the  bunch  is  positive  implying  that 
energy  is  being  extracted  from  the  optical  wave  decreasing  the  gain.  The  strong 
fields  cause  the  electron  bunch  to  move  too  far  along  the  phase-axis  in  phase- 
space  and  cause  the  electron  phase  to  shift  to  a  relative  phase  for  absorption. 
The  electrons  become  trapped  in  closed  phase-space  orbits  as  some  electrons 
can  overtake,  or  fall  behind  other  electrons  in  the  beam. 

An  estimate  of  how  much  energy  can  be  lost  by  the  electron  bunch  in  strong 
fields  can  be  made  by  looking  at  the  phase-space  diagram.  Figure  7  shows  the 
FEL  phase-space  evolution  for  weak  fields  and  low  current  where  the  pendulum 
and  optical  equations  are  solved  numerically  in  an  optical  field  a^=2  and  j=1. 
Electron  phase  velocity  v  verses  electron  phase  C  describes  the  evolution  of 
electrons  within  an  optical  wavelength.  Shown  in  the  figure  are  20  electrons;  each 
electron’s  initial  conditions  determine  the  evolution  of  its  velocity  v(t)  and  position 
C(t).  Each  electron  is  therefore  constrained  to  follow  a  particular  path  in  phase- 
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space.  The  final  phase-space  positions  of  the  electrons  are  drawn  in  (C.v).  The 
separatrix  path  separates  open  orbits  from  close  orbits.  The  height  of  the  closed- 
orbit  region  is  determined  by  the  optical  field  strength  and  ultimately  determines 
the  laser  gain  process.  The  two  plots  to  the  right  in  the  figure  represent  the 
evolution  of  the  gain  G  and  optical  phase  <p  from  the  beginning  to  the  end  of  the 
undulator. 


Figure  7.  Phase  Space  Diagram  for  Weak  Fields  with  Lx>w  Current 

The  electrons  are  initially  positioned  such  that  work  is  done  on  them;  they 
gain  energy  and  move  ahead  of  the  average  flow  at  the  top  of  the  phase-space 
picture.  Other  electrons  lose  energy  to  the  radiation  field,  and  move  back  behind 
the  average  flow  to  cause  bunching.  At  x  =1 ,  there  is  no  gain  or  phase  shift  from 
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the  uniform  spread  electron  beam,  but  as  bunching  develops  the  gain  and  optical 
field  increase.  The  maximum  energy  of  the  electron  bunch  occurs  when  electrons 
are  positioned  at  the  highest  point  on  the  separatix-axis,  and  similarly,  the 
minimum  energy  occurs  when  the  electrons  are  located  at  the  lowest  point  of  the 
v-axis.  The  total  energy  that  the  electron  bunch  can  transfer  to  the  optical  wave 
is  thereto  e  estimated  by  the  height  of  the  separatix  path,  4|ap'2  from  peak-to- 
peak. 

Figure  8  shows  the  compact  FEL  phase-space  evolution  for  weak  fields  with 
moderate  current  where  the  pendulum  and  optical  equations  are  solved 
numerically  in  an  optical  field  3^=2  and  J=5.  The  electrons  start  outside  the  closed 
orbits  and  evolve  into  a  bunch  at  C=it.  The  optical  amplitude  is  driven  by  this 
bunch  and  drives  the  gain  by  a  factor  of  100  higher  than  the  previous  case  j=1. 

vVhen  |a|=K  at  saturation,  the  change  in  the  phase  velocity  is 

Av(z)  =  —  Al  and  corresponds  to  a  change  in  the  bunch’s 

Y  2N 

energy  by  Av  =  AnNAy  I  y  =  2ti.  A  change  in  phase  velocity  by  2n  can  cause  the 
electron  bunch  to  shift  from  peak  gain  to  peak  absorption. 

Figure  9  shows  the  compact  FEL  phase-space  evolution  for  strong  field 
a^=20,  with  moderate  current  j-5.  Due  to  the  strong  fields  and  the  rapid  rate  of 
bunching  near  the  beginning  of  the  undulator,  the  gain  and  phase  begin  to  grow 
rapidly.  At  about  t=1/2,  the  optimum  relative  phase  C+4>=rt,  for  gain  occurs 
resulting  in  the  best  spatial  overlapping  of  the  electrons.  However  absorption 
occurs  as  the  bunch  evolves  further  past  this  optimum  phase  to  C+4)=0.  The  gain 
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***  FEL  Phase  Space  Evolution  *** 

j=5  a  V  =2.6  N=49 
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Figure  8.  Phase>Space  Diagram  for  Weak  Fieids  with  Moderate  Current 

is  seen  to  first  peak,  and  then  decrease  as  the  overbunched  electron  beam  drives 
the  optical  wave  and  continues  to  evolve  so  that  energy  flows  back  into  the 
electron  beam.  The  final  gain,  G=10%,  is  significantly  reduced  from  the 
theoretically  possible,  G=0.135j,  in  weak  optical  fields. 

As  anticipated,  the  FEL  interaction  in  strong  optical  fields  demonstrates 
significantly  different  characteristics  upon  saturation.  The  gain  spectrum,  G(vo). 
changes  shape  in  the  presence  of  strong  fields.  As  discussed  in  the  gain 
mechanism  section,  the  gain  spectrum  in  weak  fieids  and  low-current  is  anti¬ 
symmetric  about  the  resonance  wavelength  with  a  peak  gain  G=0. 135j,  at  phase 
velocity  Vo=2.6  and  has  a  gain  bandwidth  of  roughly  Avo=ji.  Similarly  the  peak 
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***  FEL  Phase  Space  Evolution  *** 
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Figure  9.  PhaseoSpace  Diagram  for  Strong  Fields  with  Moderate  Current 

absorption  occurs  at  the  phase  velocity  Vo=-2.6  with  a  value  of  G=-0.135j.  For 
j=5,  the  gain  spectrum  is  shifted  slightly  to  the  left  away  from  resonance.  As  the 
initial  field  strength  increases,  the  peak  gain  decreases  and  the  wavelength  at 
which  peak  gain  occurs  also  increases  away  from  resonance.  The  gain  spectrum 
width  becomes  broader  and  distorts  as  the  FEL  saturates.  The  gain  is  further 
reduced  as  the  bunched  electrons  move  out  of  resonance  and  lose  energy  to  the 
optical  field.  These  effects  are  shown  in  Figure  10  where  the  gain  spectrum  G(vJ 
for  moderate  current  y=5  and  initial  field  strength  8^=0  is  increased  to  30.  The 
phase  velocity  that  gives  peak  gain  increases  away  from  resonance  as  the  field 
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strength  inaeases.  The  gain  describing  the  interaction  after  each  pass  in  the 
FEL  oscillator  causes  the  optical  field  to  grow  a  little  stronger  as  can  be  seen 
along  the  a^-axis.  As  the  power  grows  over  many  passes,  typically  several 
hundred  passes  are  required  to  evolve  from  the  weak-field  to  the  strong-field 
regime  near  saturation,  the  optical  wavelength  is  observed  to  shift  to  longer 
wavelengths  corresponding  to  a  larger  Initial  phase  velocity. 


In  (l-i-G) 


Figure  10.  The  Gain  Spectrum  G  fvj  for  Moderate  Current  with 
Increasing  Optical  Field  Strength  a^. 


E.  PRODUCTION  OF  HIGHER  FREQUENCY  HARMONICS  IN  FELs 

The  frequency  of  radiation  emitted  by  a  relativistic  electron  traveling  through 
a  linearly-polarized  undulator  is  given  by  nw  » <i),{  1  )  /  2y*.  Here,  n=1 ,3,5,.., 
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is  known  as  the  harmonic  number.  The  fundamental  is  defined  as  the  frequency 
produced  when  n=1  [Ref.  11].  We  address  an  overview  of  the  harmonic 
mechanism  to  provide  an  intuitive  understanding  of  their  origin  so  specific  compact 
designs  that  operate  on  harmonics  can  be  explored  in  Chapter  IV.  Using 
numerical  methods.  Chapter  IV  will  also  examine  the  advantages  of  using 
harmonics  in  compact  FELs.  Here,  dynamics  of  relativistic  electrons  in  the 
presence  of  light,  stimulated  emission,  are  generalized  to  describe  higher 
harmonics  of  the  optical  field  [Ref.  6.12].  Harmonics  are  produced  by  extra 
oscillations  of  the  electron  that  occur  at  each  undulator  period.  The  energy 
transferred  between  the  electron  beam  and  light,  given  in  (10).  is  extended  to 
include  this  motion. 

Harmonics  result  from  a  coupling  between  electrons  and  light  due  to  an 
alteration  of  the  energy  transfer  mechanism.  In  addition  to  the  transverse 
sinusoidal  motion  of  electron  trajectories  in  the  linearly-polarized  undulator,  a  small 
periodic  longitudinal  motion  exist  that  cause  spontaneous  emission  and  gain  in 
odd  multiples  of  the  fundamental,  n(o=nkc.  The  harmonic  coupling  is  due  to  the 
time  electrons  spend  in  periodic,  longitudinal  motion  while  transferring  energy  to 
the  optical  wave.  A  weighted  coupling  between  electrons  and  light  is  expressed 
by  Bessel  functions  which  modify  the  undulator  parameter  K.  The  modified 
undulator  parameter  provides  a  convenient  method  of  incorporating  harmonics  in 
the  optical  wave  equation,  electron  pendulum  equation  and  gain  analysis.  An 
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outline  of  how  this  modified  undulator  parameter  is  derived  and  incorporated  into 
(10)  is  given  below. 

The  ideal  undulator  field  is  given  by 

5,  =  S(0.sin(  V).0)- 

The  components  of  the  Lorentz  force  equations  including  this  new  undulator  field 
are 

a  <  ’’I’  > '  ^  ^  '  7^(  “■ 

^  =  0  =►  Y  =  constant .  (45) 

The  X  and  /-velocities  are  found  by  integrating  the  respective  component  equation: 

p,  =  )!^cosk^z  *  6  p^ 

P^  =  constant  +  ftp,  (47) 

where  K=  eBk.l2nmc^  and  ftp,  and  ftp,  are  constants  of  integration.  Using 
P^  =  1  -  Y'^  -  Px  -  Px  .(46)  and  (47),  and  assuming  perfectly  injected  electrons  with 
ftp,  =  ftpy  =  0,  the  electron  z  motion  is  given  by 


36 


Pf  =  1  -Y'^[(1  *  K^cos{2k^z)  ' . 


(48) 


At  z=0,  the  injected  z  velocity  is 

p?  u*p;(o)=i  -ymi  *2«’) 


PJ,..*P.(0)  =  Y 

for  perfect  injection. 

Now  (48)  can  be  written 

P'  =  Pz{0)-j^j"-(^|'cos(2V) 


(49) 

(50) 


(51) 


and  perturbation  expansion  in  (Klyf  is  now  possible.  Defining 
Pr  ( 0  *  P,  ♦  A  p,  ( f ) .  the  2  motion  to  order  (Klyf  is  derived  from  integrating  z  to 
lowest-order  to  give: 


Pf  =  P*-2j5|^cos(2/c„z). 


(52) 


To  lowest  order  in  (Klyf.  [P,  =  2/c=p,  =  1  -1  '2(1  ^  1^)]  so  that  the  electron 
position  can  be  written  as  the  expansion  +  P,cf  + The  integrated  z 
motion,  becomes 


n  -  ( ^  ^ f  sin(  ( 2k,z, )  -  ) . 


(53) 


The  first  two  terms  express  the  slowly  evolving  motion  of  the  electron  along  the 
undulator-axis.  The  last  term  expresses  the  fast,  periodic,  longitudinal  motion  on 
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the  optical  wavelength  scale  and  causes  higher  harmonics.  This  oscillatory  motion 
occurs  once  per  undulator  period.  Variations  of  the  fundamental  electric  field 
observed  by  the  electron  as  it  travels  through  a  undulator  period  causes  the 
nonuniform  axial  motion.  The  fast  and  slow  motions  can  now  be  expressed  as 

z(  f)  =  I  ^j^sin(2o>^f)  where  z^  =  0  and  =  1  • 

Defining  $  =  /  2  ( 1  -  ) .  (54)  can  be  expressed  as 

z(0  =  3,cf-|sin(2a),0.  (55) 


Substitution  of  (55)  into  (1 0)  and  recalling  C  =  (k  k^}z  -  oit  is  the  electron  phase, 
the  electron  energy  equation  is  written  explicitly  in  terms  of  fast  and  slow  variables: 

Y  =  ( cos(  nc  -  (p  -  ( n  -  1  )o),f  -  r75sin(  2ioJ) ) 

*  cos(  nC  *  <p  -  ( n  -  1  )<jJ  -  r7$sin(  2(oJ)  ] 

Using  the  generating  function  to  expand  sinusoidal  terms  in  cylindrical  Bessel 
functions  and  averaging  over  one  undulator  period,  (56)  is  written  as 


Y  = 


eBtJ^EL 
2y  mcf 


cos  ( nC 


<P ) 


where 


=  ^( -1  ,2{ no  - „2(nU]  for  n  =  1,3,5...  . 

The  product  of  the  undulator  constant  and  the  Bessel  functions  provides  a 
convenient  mathematical  notation  to  compare  FEL  operation  in  any  selected 
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harmonic.  By  using  the  theory  already  developed  in  the  preceding  sections  the 
generation  of  harmonics  can  be  expressed  by  simply  replacing  /<-  $  )• 
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III.  CONCEPT  OF  A  COMPACT  FEL 


To  increase  the  accessibility  of  FEL’s  to  a  broad  range  of  users  a  reduction 
in  size  and  complexity  of  FEL  systems  must  be  pursed.  The  largest  component 
of  a  FEL  system  is  usually  the  electron  accelerator  due  to  the  large  electron 
energy  required.  The  relatively  large  electron  beam  energy  along  with  high 
current  requires  significant  radiation  shielding  for  personnel  safety  and  therefore 
increasing  the  overall  weight  of  the  FEL  system  [Ref.  8].  To  reduce  these 
requirements,  FELs  that  operate  on  harmonics  have  been  proposed  to  allow 
using  lower  energy  electron  beam.  This  effect  can  be  examined  by  considering 
a  fixed  undulator  design,  and  referring  to  the  resonance  condition  given  in  [Ref. 
12],  ^X  =  X^{^  +/<^)/2y^.  where  the  harmonic  number  n  is  included.  The 
undulator  design  depends  on  the  period  and  undulator  parameter  K  and 
therefore  assumed  to  be  constant.  For  a  desired  optical  wavelength,  the  only 
variables  in  the  resonance  equation  are  the  Lorentz  factor  y  and  harmonic  number 
n.  By  increasing  n,  we  see  y  can  be  reduced  accordingly,  provided  K  and  k 
remain  fixed.  Since  y  is  proportional  to  the  electron  beam  energy,  using 
harmonics  can  provide  a  viable  method  of  designing  compact  FELs  which  require 
lower  energy. 
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A.  DEVELOPMENT  OF  THE  PHOTO-CATHODE  INJECTOR 

Preliminary  design  studies  have  been  conducted  at  Los  Alomos  National 
Laboratory  (LANL)  on  an  infrared  compact  RF-linac  driven  FEL  [Ref.  8]  that  relies 
upon  two  technological  advances: 

•  a  very  bright  electron  beam  from  a  photo-cathode  injector,  and  a 

•  pulsed-wired  technique  for  short-period  high-field  undulators. 

The  photo-cathode  injector,  designed  by  Sheffield  and  Fraser  at  LANL,  is  a 
compact  device  used  to  generate  the  FEL’s  electron  beam  [Ref.  13,14],  The 
device  uses  a  light-activated,  photo-emissive  electron  source  in  order  to  provide 
excellent  control  of  the  electron  distributions.  Electron  distribution  functions 
determine  the  quality  of  the  electron  beam,  and  play  a  significant  role  in  the  small- 
signal  gain  that  is  available  from  a  given  undulator.  Beam  quality  can  be 
described  by  two  effects,  emittance  and  energy  spread  [Ref.  5].  Both  effects 
reduce  FEL  performance  by  degrading  the  bunching  process  described  earlier. 
The  transverse  emittance  is  used  to  describe  the  imperfect  FEL  injection  of 
electrons  off-axis  in  either  position  or  angle.  Since  the  transverse  coordinates  of 
the  electrons  in  the  FEL  are  (x,y).  we  must  consider  electron  injection  off-axis  to 
include  both  dimensions.  Figure  1 1 .  however,  shows  the  electron’s  injection  only 
in  the  /-direction.  A  similar  schematic  holds  for  the  x-direction. 

An  imperfectly  injected  electron  has  a  decrease  in  phase  velocity,  and 
random  spreads  in  {  /^ ,  0^ )  and  ( x^ ,  0,  ) .  These  effects  will  begin  to  decrease 
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Figure  11.  Schematic  of  Injection  of  imperfect  Electrons  into  Undulator. 
Injection  of  imperfect  off*axis  into  the  undulator  due  to  spread  of  position 
and  angle  6,^. 

the  PEL'S  ability  to  bunch  electrons  causing  a  decrease  in  gain.  Assuming  the 
motion  of  the  electrons  in  the  two  x  and  y  transverse  directions  is  independent,  a 
transverse  phase-space  plot  can  be  constructed.  Taking  the  electron  distribution 
to  be  a  product  of  uncorrelated  Gaussians  in  angular  spreads  and  positions,  the 
transverse  emittance,  and  ^  are  defined  by  the  corresponding  phase-space 
areas.  See  Figure  12. 

(57) 

An  ellipse  drawn  around  the  electron  particle  distribution  defines  an  area  in 
y-B  phase  space  that  characterizes  the  transverse  emittance.  Similar  results  hold 
in  x-d  phase  space. 
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Figure  12.  Electron  Particle  Phase>Space  Distribution. 

An  ellipse  drawn  around  the  electron  particle  distribution  defines  an  area  in 
y-6^  phase  space  that  characterizes  the  transverse  emittance.  A  Gaussian 
distribution  of  positions  occur  in  y  with  an  uncorrelated  Gaussian 
distribution  of  angles  in  6^.  Similar  results  hold  in  x-d,  phase  space. 


The  emrttance  growth  of  the  beam  is  reduc  1  by  using  the  photo-injector 
design,  because  the  electron  bunch  produced  fr^m  the  cathode  is  rapidly 
accelerated  to  relativistic  energies  in  a  single  cavity,  thereby  reducing  significantly 
the  electron  beam  transport  at  low  energies  [Ref.  13,14].  The  desired  low 
emittance  is  also  required  for  the  electron  beam  to  propagate  through  the  small 
gap  in  the  undulator,  and  the  optical  mode  in  the  resonator  cavity.  Brightness,  a 
property  of  the  electron  beam,  is  proportional  to  the  current  and  inversely 
proportional  to  the  beam  quality  [Ref.  1 5].  This  implies  emittance  must  be  reduced 
in  order  to  obtain  an  intense,  bright  electron  beam. 
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B.  TECHNOLOGICAL  ADVANCES  IN  MICRO-UNDULATOR  DESIGNS 


Warren  and  Feldman  [Ref.  16]  discuss  the  design  and  construction  of  iron- 
free  electromagnets  driven  by  pulsed  high-currents.  The  micro-undulator  can 
generate  the  high  magnetic  fields  needed  to  reduce  the  period,  and  still  maintain 
optical  gain  at  reasonable  values.  These  devices  have  periods  of  less  than  one 
centimeter  and  can  be  used  to  generate  light  at  short  wavelengths  with  electron 
beams  of  only  moderate  energy.  White  reducing  the  period,  the  designer  requires 
the  undulator  parameter,  K,  to  be  less  than  one  in  order  to  extract  the  maximum 
gain  from  the  FEL.  Pulsed  currents  in  copper-like  conductors,  that  produce 
current  densities  in  excess  of  about  1 0^  A/cm^,  can  be  utilized  to  accomplish  this 
requirement.  Limits  of  this  technology  involve  undulator  construction  errors, 
material  stresses  and  electron  beam  alignment  issues  [Ref.  17].  Compared  to 
other  undulator  technologies,  this  design  has  the  greatest  potential  to  reach  short 
FEL  wavelengths  with  reduced  electron  beam  energy. 

The  harmonic  mechanism  in  FELs,  discussed  in  the  last  chapter,  is  the  result 
of  the  fast,  periodic  axial  motion  of  the  electrons.  The  effects  of  harmonics  can  be 
m.anipulated  by  designing  micro-undulators  that  operate  with  periods  which  are 
multiples  of  the  fundamental  period.  Re-writing  the  resonance  condition  for 
harmonics  as  X  =  A^(  1  +  K^)  12  y^n,  it  can  be  shown  harmonics  of  the  undulator’s 
period  can  be  utilized  to  design  micro-undulators.  These  designs  would  have 
fewer  periods  and  still  have  the  requirements  mentioned  earlier. 
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The  top  schematic  of  Figure  13  illustrates  an  FEL  operating  on  the 
fundamental,  n  =  1,  for  a  fixed  undulator  and  fixed  electron  beam  energy.  In  the 
fundamental,  the  undulator  has  a  period  undulator  parameter  K,  and  initial 
electron  velocity  corresponding  to  a  Lorentz  factor,  y.  The  fundamental  optical 
wavelength,  k,  is  found  from  the  harmonic  resonance  condition  with  n  =  1.  The 
middle  schematic  illustrates  the  same  undulator  design  and  electron  beam  energy, 
but  the  third  harmonic,  n=3,  is  used.  The  resulting  wavelength  is  kl3  or  one  third 
the  fundame.ntal  wavelength  satisfying  the  harmonic  resonance  condition  with  n=3. 
This  is  the  conventional  use  of  harmonics  to  reduce  the  optical  wavelength.  A  new 
method  of  utilizing  harmonics  proposed  by  R.  Warren,  is  shown  at  the  bottom. 
Using  the  same  electron  beam  energy  and  undulator  parameter  K,  the  optical 
wavelength,  k,  is  acquired  from  the  harmonic  resonance  condition  when  the 
undulator  period  is  increased  to  three  times  the  fundamental’s  undulator 
wavelength  3  k^.  The  advantage  of  this  design  is  that  the  micro-undulator  FEL  can 
reach  short  optical  wavelengths  with  a  reasonable  undulator  period  and  a  low 
energy  electron  beam.  The  next  chapter  discusses  in  detail  three  specific  micro- 
undulator  designs  that  utilize  harmonics  in  this  fashion. 
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Figure  13.  Schematics  of  Compact  FEL  Operating  on  Optical  and 

Undulator  Harmonics. 
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IV.  EVALUATION  OF  COMPACT  FEL  OPERATION 


A.  HIGHER  HARMONIC  GENERATION  IN  COMPACT  FELs 

There  is  currently  interest  in  lasing  on  higher  harmonics  of  the  fundamental 
frequency  in  order  to  reach  optical  wavelengths  in  the  range  of  0.4  to  0.5  microns 
with  an  electron  beam  energy  of  about  20  Mev.  The  use  of  harmonics  at  short 
wavelengths  allows  a  realistic  undulator  parameter,  since  the  undulator’s  length 
is  shorter,  the  number  of  periods  is  less  and  longer  periods  would  be  used.  The 
consequence  would  be  a  micro-undulator  design  that  would  be  easier  to 
fabricate.  Shorter  undulators  would  also  minimize  the  penalty  from  emittance  and 
energy  spread  inherent  in  all  realistic  electron  beams  [Ref.  18].  For  a  fixed  optical 
wavelength,  A,  and  an  undulator  design  using  higher  harmonics,  the  energy  of 
the  electron  beam,  ymc^,  can  be  reduced  thereby  lowering  the  radiation  shielding 
requirements  for  personnel.  Operation  over  a  wider  range  of  wavelengths,  and  the 
increase  in  optical  power  at  saturation  provide  additional  benefits  of  using  higher 
harmonics  in  compact  FELs  [Ref.  12],  The  gain  that  can  be  achieved  using 
various  harmonics  depends  on  electron  beam  properties,  beam  quality  factors, 
and  the  specific  micro-undulator  design  chosen.  To  evaluate  gain  on  harmonics 
for  new  FEL  design  systems,  it  is  necessary  to  utilize  computer  codes  that 
represent  the  mathematical  formulation  of  the  FEL  equations  discussed  in  Chapter 
II. 


47 


B.  FEL  SIMULATION  METHODS 


Various  numerical  simulation  methods,  essential  to  relating  the  three  major 
technologies  of  optics,  relativistic  electron  beams  and  undulator  magnets  have 
been  developed  to  analyze  and  design  RF-linac-driven  FELs  [Ref.  19].  Numerical 
calculations  using  five  methods  that  model  FEL  laser  physics  with  different  levels 
of  sophistication  and  mathematical  formalism  are  conducted  on  three  micro- 
undulator  designs.  The  designs  are  explored  to  demonstrate  the  use  of  higher 
harmonics  in  order  to  reach  optical  wavelengths  with  lower  electron  beam  energy. 
These  preliminary  compact  FEL  designs  were  proposed  by  R.  Warren  at  LANL 
The  design  method  is  based  on  Dattoli’s  parameterized  gain  formula  [Ref.  20]  and 
approximate  solutions  for  the  magnetic  fields  of  single  and  double-current  carrying 
helical  windings  [Ref.  21].  The  intent  of  the  design  code  is  to  accurately  model 
gain  characteristics,  and  provide  a  basis  within  which  micro-undulator  parameters 
can  be  developed  in  order  to  optimize  compact  FEL  operation.  The  assumptions, 
approximations,  regions  of  validity  and  definitions  of  important  parameters  for  the 
methods  used  in  the  analysis  will  be  discussed  through  this  chapter. 

C.  COMPACT  FEL  DESIGN  PARAMETERS 

The  electron  beam  current  and  energy  spread,  and  the  optical  wavelength 
are  the  same  for  all  three  micro-undulator  designs  evaluated.  The  electron  beam 
has  current  /=400  A,  and  Lorentz  factor  y=40  corresponding  to  an  energy  of 
about  20  MeV.  The  beam  emittance  is  varied  from  g=A./4  to  e=7k,  while  the 
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Gaussian  energy  spread  varies  from  0%  to  0.8%.  The  optical  wavelength  is  A =0.4 
microns  with  the  optimum  Rayleigh  range  =  LI  /fl”  where  L  is  the  length  of  the 
undulator.  Two  of  the  three  undulators  are  designed  to  be  operated  on  higher 
harmonics  with  undulator  periods  in  the  range  of  X^=0.3  to  0.5  cm.  The  values 
of  the  harmonic  number,  number  of  periods,  period,  and  undulator  parameter  are 
given  for  the  three  micro-undulators  in  Table  I. 


TABLE  I. 


FUNDAMENTAL 

THIRD  HARMONIC 

NINTH  HARMONIC 

n 

1 

3 

9 

N 

413 

49 

19 

A,(cm) 

.013 

0.33 

0.54 

K 

0.018 

0.408 

0.540 

D.  GAIN  DEGRADATION  DUE  TO  BEAM  QUAUTY  AND  THE  FILUNG 
FACTOR 

An  analysis  of  beam  quality  and  the  filling  factor  ^  is  necessary  for  a 
discussion  of  the  FEL  gain  evaluation.  Emittance  and  energy  spread  effects  that 
describe  beam  quality,  introduced  in  Chapter  III,  are  discussed  in  greater  detail  to 
show  important  FEL  design  limits.  The  filling  factor  9^,  defined  as  "the  ratio  of  the 
electron  beam  size  nr/"  to  the  "optical  mode  size  nw/,"  is  used  to  define  a  limit 
on  the  electron  beam  area.  [Ref.  5] 
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Gain  in  the  FEL  depends  on  details  of  the  phase  velocity  distribution  where 
the  ability  to  bunch  electrons  is  diminished  in  a  different  way  for  each  shape  [Ref. 
22,23].  The  exact  distribution  function  of  the  electrons  transverse  motion  depends 
on  the  construction  and  operation  of  the  accelerator  and  usually  is  not  known 
exactly.  The  distribution  for  initial  electron  energies  perfectly  injected  into  the 
undulator  is  represented  by  the  normal  distribution  function  [Ref.  23] 


/g(v) 


ei<p(-(v-v,)^/24) 


for  all  V 


(58) 


where  fgfvj  =  Gaussian  distribution 

V  =  phase  velocity 

V  =  peak  phase  velocity 

Oq  =  4nN  A  y/y  is  the  standard  deviation  of  v  away  from  the  peak 
phase  velocity  v^. 


A  Gaussian  distribution  of  angles  about  the  z-axis  is  given  by  the  exponential 


distribution  function  [Ref.  23] 


/0(V) 


exp{-(v^-v)/  09] 
00 


for  v<v^ 


/( V  )  =  0  for  \>\^ 


(59) 


where  =  4nNY*(  A0/1  ■►^*) 

A0  =  standard  deviation  of  angular  disstribution 
Vo  =  phase  velocity  for  electrons  entering  on-axis. 

It  is  important  to  note  that  for  any  angular-spread  of  phase  velocity,  the 

longitudinal  velocities  of  these  electrons  are  decreased.  The  phase  velocity  is 

skewed  only  in  the  negative  direction  indicating  that  the  electron’s  energy  can  only 

be  decreased  regardless  of  the  initial  angle  relative  to  the  z-axis. 
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The  decrease  in  phase  velocity  for  an  imperfectly  injected  electron  is  given 
by  [Ref.  5], 

(60) 

where  Av  =  change  in  electron’s  phase  velocity,  and 
0^^  =  electron’s  initial  injection  angle. 

Any  change  in  the  electron’s  average  z-velocity  causes  a  change  in  the  electron’s 
phase  velocity  as  given  by  the  above  expression.  When  the  electrons  are  injected 
off-axis  in  either  position  or  angle,  they  experience  a  field  strength  that  is  stronger 
than  the  field  strength  on  the  central-axis.  This  effect  causes  the  electrons  to  be 
focused  back  toward  the  center  line,  and  execute  a  number  of  betatron 
oscillations,  =  NK/  y ,  in  the  transverse  x-y  plane.  The  betatron  oscillations  result 
from  transverse  excursions  of  electrons  that  have  a  constant  of  motion  defined  by 
a  harmonic  potential  [Ref.  5].  In  the  situation  where  the  electron  beam  is  guided 
into  this  betatron  motion,  it  is  desirable  to  match  the  angular  and  position  spreads 
so  that  the  beam  does  not  focus  or  expand  as  it  travels  along  the  z-axis.  A  smooth 
propagation  of  the  electron  beam  with  minimum  spread  in  the  phase  velocities  is 
achieved  by  requiring  that  the  angular  and  position  spreads  be  matched  [Ref.  5]; 

Kk,r,  •  y0,  (61) 

where  0^  is  now  a  spread  in  angles. 

The  spread  of  electron  phases  accumulated  over  the  undulator  length  in  a 


matched  beam  is  given  by  [Ref.  5] 


(62) 


Av=J^K^k!r!. 

1  >  /<2 

If  we  require  the  phase  spread  not  exceed  an  optical  wavelength  at  the  end  of  the 
undulator,  AzlA^I  the  limit  for  is  found. 

Relating  Av  to  Azik  at  the  end  of  the  undulator  we  first  find  A(i;  =  Av  at  t  =  1 . 
since  V  =  C-  Recalling  A(,=kAz,  the  phase  spread  can  be  written  as  Av=kAz. 
Eliminating  Az  from  Azlk>^  and  A\)=kAz,  and  using  (62)  we  find  the  phase  spread 
limit  can  be  written  as 


AnN  ^  ^  (63) 

Using  the  resonance  condition,  =  NK/y  and  L  =  Nk^,  the  limit  for  is  found 
to  be 


nLX 

An^Nf 


(64) 


This  resulting  limit  on  the  electron  beam  area  is  to  insure  that  the  phase  spread 
does  not  degrade  beam  quality. 

An  important  limit  of  the  electron  beam  radius  constrains  the  electron  beam 
to  be  inside  the  optical  mode  w^.  The  fraction  of  the  beam  outside  the  mode 
cannot  participate  in  the  gain  process.  The  expansion  of  a  coherent  beam  with  a 
Gaussian  intensity  profile  in  the  fundamental  mode  is  given  by  [Ref.  24] 
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W^{z)  =  kV< 


2 

1  + 

\x{z-L/2) 

nWo 

\  } 

(65) 


where  w(z)  describes  the  beam  contour  at  position  z,  and  the  waist  spot  size 
at  z=0. 


The  Rayleigh  range,  defined  as  =  n  A,  is  the  distance  where  the  optical 
transverse  mode  expands  to  twice  the  size  compared  to  waist  size  in  the  undulator 
[Ref.  5].  See  Figure  14.  The  Rayleigh  range  is  optimized  by  minimizing  the  optical 
mode  volume  around  the  electron  beam.  The  optical  mode  volume  is  calculated 
by  integrating  the  beam  intensity  profile  along  the  entire  length  of  the  undulator 
[5], 


i/2 


V^2  \  n 


w! 


i.£! 

z! 


dz  =  XL 


O  / 


(66) 


where  X  I z^. 

The  two  limiting  cases,  z^->  0  and  z„— <»,  each  result  in  an  optical  mode 
volume  that  tends  to  infinity.  In  the  first  case,  a  small  waist  is  pinched  in  the  middle 
of  the  optical  mode,  and  rapidly  expands  due  to  diffraction.  An  infinite  volume 
means  that  light  at  the  ends  of  the  undulator  would  be  far  from  the  electron  beam 
and  would  not  be  significantly  amplified.  In  the  second  case,  a  plane  wave  would 
develop  where  all  of  the  optical  mode  is  large.  The  optical  mode  would  be  far 
away  from  the  electron  beam  at  the  end  of  the  undulator  and  at  the  waist,  as  well. 
Both  these  cases  are  unsuitable  for  any  physical  solution  therefore  their  must  be 
a  optical  mode  volume  minimum  between  0  and  ». 
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Figure  14.  Schematic  for  Diffraction  of  Light  Wave  In  Undulator. 
Schematic  for  a  Gaussian  beam  in  lowest-order  expanding  due  to  diffraction 
spreading  as  it  propagates  away  from  the  waist  region  located  at  center  of 
the  undulator. 

We  find  the  value  of  that  make  the  volume  of  the  optical  mode  a  minimum 
(Ref.  Setting  the  derivative  of  (66)  to  0 


the  optimum  Rayleigh  range  is  found  to  be  -  IZ/ii*.  By  making  the  Rayleigh 
range  equal  to  about  a  third  to  a  fourth  the  undulator.  the  minimum  optical  mode 
volume  is  found.  This  characteristic  volume  provides  the  optical  mode  for 
maximum  gaia 

For  the  optimum  Rayleigh  range,  s  gives  a  limit  to  the  electron  beam 

area 
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Lk 


(68) 


Tcr|  ^ - 

When  this  condition  is  satisfied,  the  small-signal  gain  is  independent  of  the  electron 
beam  size,  and  depends  only  on  the  current  contained  within  the  optical  mode. 
Emittance,  the  product  of  the  rms  values  of  the  angular  and  radial  position 
spreads,  provides  the  fundamental  limit  to  beam  size  and  determine  a 
corresponding  spread  in  the  longitudinal  electron  phases.  Comparing  the 
restrictions  (62)  and  (66)  shows  that  for  ,  the  only  effect  of  emittance  occurs 
when  the  electron  beam  is  larger  than  the  optical  mode,  or  simply  F  .  This 
restriction  simplifies  the  analysis  of  any  PEL 

E.  METHODS  OF  ANALYSIS 

1.  Three-Dimensional  Simulation  Codes 
a.  FELEX 

FELEX,  further  abbreviated  to  FEX,  numerically  calculates  three 
dimensional  particle  motion  and  optical  diffraction  effects  [Ref.  25].  The  computer 
code  is  based  upon  the  Monte  Carlo  technique  that  simulates  three  dimensional 
FEL  physics  by  following  orbits  of  simulation  electrons;  2,560  are  used  in  our 
analysis  to  determine  the  optical  field  interaction.  Numerical  solutions  to  the  FEL 
equations  (5a),  (5b),  and  (6)  are  found  using  three  major  algorithms  [Ref.  25],  (i) 
a  set  of  ordinary  differential  equations  describing  the  electron  motion,  (ii)  a 
parabolic  partial  differential  equation  governing  the  optical  field,  and  (iii)  the 
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connection  between  the  optical  field  and  the  electron  motion  through  a  source 
term  that  is  a  function  of  the  electron  current  density.  The  currents  from  the  set 
of  electron  equations  provide  the  source  term  in  Maxwell’s  equations  and  connect 
the  electron  motion  to  the  optical  field  [Ref.  25].  Averaging  the  relativistic  Lorentz 
force  equation  (1)  over  the  spatial  period  of  the  undulator,  the  numerical  solution 
is  acquired  that  incorporate  in  the  dynamics  (i)  finite  emittance,  (ii)  energy  spread, 
and  (iii)  linear  undulator  focussing  of  the  electron  beam  in  the  undulator  field.  The 
program  was  set  up  for  curved-faced  pole  undulators  where  equal  focusing  in  the 
X  and  y-directions  is  obtained.  This  effect  causes  the  electron  beam  to  remain 
round  as  it  travels  through  the  undulator  instead  of  being  focused  only  in  one 
direction,  and  thereby,  developing  an  elliptical  shape.  In  short  undulators,  which 
are  considered  here,  gain  calculations  are  not  substantially  affected  by  the  fact  that 
the  simulation  code  uses  curved-faced  poles  instead  of  flat-faced  poles.  These  are 
single  wavefront  calculations  that  do  not  consider  the  pulse  nature  of  both  the 
electron  and  optical  beam.  Optical  gain  and  refractive  effects  are  "self-consistently" 
calculated  as  the  light  interacts  with  electrons  inside  the  undulator.  A  finite  width 
of  the  optical  wavefront  near  the  center  of  the  optical  pulse  length  is  modeled.  The 
transverse  profile  of  the  electric  field  is  set  equal  to  a  normal  Gaussian  mode 
characterized  by  a  wavelength,  focal  position,  and  a  Rayleigh  range.  The  optical 
field  at  the  exit  of  the  undulator  is  propagated  through  the  optical  resonator  to  the 
entrance  of  the  undulator  where  the  next  electron  pulse  interacts  to  simulate  an 
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FEL  oscillator.  The  evolution  of  the  optical  field  is  found  by  solving  the  paraxial 
wave  equation  [Ref.  25,26],  since  FEX  models  the  optical  components  of  the 
resonator  by  thin  lens.  The  paraxial  wave  equation,  derived  from  the  formulation 
for  the  optical  wave  equation  discussed  in  Chapter  II,  has  an  inhomogeneous 
driving  term  that  is  a  result  of  the  stimulated  emission  of  the  electron  beam  when 
placed  in  the  optical  field.  An  extension  of  the  FEX  code  was  made  at  LANL  to 
incorporate  a  model  of  coherent  spontaneous  emission  of  harmonics  in  the  FEL 
oscillator  operating  at  high  power  at  the  fundamental  wavelength  with  a  plane 
polarized  undulator.  An  additional  modification  allows  calculations  of  undulator 
field  errors  due  to  manufacturing  flaws  in  magnetic  material  or  imperfect  alignment 
during  construction  [Ref.  26]. 

b.  Wavefront  simulation 

A  wavefront  simulation,  WFS,  numerically  solves  the  3-D  parabolic 
wave  equation,  together  with  the  pendulum  equation  [Ref.  5].  This  method 
includes  most  of  the  FEL  effects  incorporated  into  FEX  except  that  a 
monoenergetic  electron  beam  is  assumed.  When  the  gain  is  sufficiently  high,  the 
electron  beam  can  distort  the  optical  mode  into  a  combination  of  resonator  modes 
[Ref.  5].  This  code  accounts  for  distortion  by  self-consistently  calculating  the 
optical  field  for  high  gain  in  3-D  without  using  the  filling  factor.  Harmonic  coupling 
and  strong  fields  are  included  in  the  code. 
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2.  One-Dimensional  Simulation  Code 

Numerical  solutions  of  the  wave  and  pendulum  equations,  SIM,  for  both 
strong  and  weak  fields  are  found  with  approximately  1000  simulation  electrons. 
This  analysis  includes  the  filling  factor,  harmonic  coupling  and  electron  beam 
degradation  effects.  The  evaluation  of  the  FEL  interaction  is  described  by  the 
dimensionless  current  j  which  determines  the  response  of  the  optical  wave  to 
bunching  in  the  beam  [Ref.  5]. 

3.  Analytical  Methods 

a.  Integral  Equation 

Combining  the  electron  pendulum  and  optical  wave  equations  in 
weak  fields,  an  integral  equation,  INT,  governs  the  evolution  of  the  optical  field 
[Ref.  5,22].  This  development  allows  the  evaluation  of  gain  degradation  from  a 
Gaussian  electron  beam  distribution  in  both  the  high  and  low-gain  regimes  where 
the  individual  electron  phase  dependence  is  eliminated.  Coupling  to  higher 
harmonics  is  also  incorporated  in  the  integral  equation  which  is  solved  numerically 
on  a  small  computer.  The  transverse  optical  mode  remains  fixed  during  this  FEL 
interaction.  However,  fixed  optical  mode  coupling  is  incorporated  by  including  the 
filling  factor. 

b.  Small-Signal  Gain  Formula 

A  modified  small-signal  gaim  formula,  DAT,  which  incorporates 
harmonics  and  the  breeding  of  the  electron  beam  due  to  emittance  and  energy 
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spread  was  also  used  in  developing  the  gain  plots.  Dattoli's  parametrized  gain 
formula  [Ref.  20]  is  given  in  cgs  units  by 


G  -  19  X  ^  Ml  [yg  ,  )2 


1  -  a:*  V 


(1  ♦K/3)(1 


(69) 


where  G  =  single-pass  gain, 

N  =  number  of  undulator  periods, 

n  =  harmonic  number, 

K  =  rms  undulator  parameter, 

I  (A)  =  peak  current  in  amperes, 

Y  =  Lorentz  factor, 

Jq  ,  =  Bessel  function  =  [  J(  /?  -  1 )  /  2  (  nO  -  J{n  +  ^ )  /2{nl)]  where 
i  =  1/2  K^n  +/C2. 

=  coupling  parameter, 

/ij,  y  =  broading  parameter  due  to  emittance, 
n,  =  broading  parameter  due  to  energy  spread. 

The  parameterization  of  electron  beam  characteristics  contained  in  (69)  has  been 

obtained  by  a  polyncmial  expansion,  and  best  fit  with  num.erical  data  [Ref.  27,28]. 

The  coupling  parameter  is  inversely  proportional  to  the  rms  longitudinal  bunch 

length  of  the  electron  pulse  and  does  not  significantly  contribute  to  degrading  the 

gain  since  for  the  cases  considered.  For  an  ideal  linear  undulator  the 

magnetic  field  is  constant  along  the  x-direction.  By  ideal  we  mean  undulator 

magnets  that  extend  to  infinity  in  the  x-direction.  In  this  case,  electrons  are  not 

focused,  but  continue  their  motion  away  from  the  undulator-axis.  For  the  ideal 

linear  undulator  the  horizontal  emittance  can  therefore  be  neglected,  since  it  does 

not  contribute  to  gain  degradation.  Dattoli’s  parametized  gain  formula  was 

modified  by  D.  Warren  to  include  these  assumptions  [private  communication].  The 

modified  small-signal  gain  formula  (DAT)  is  written  as 


59 
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G  =  19x10-‘*Ar2„2 


1  Y 


1  * 


^  Y  I 


1  *1.7(£Sr 


(70) 


where  ES  =  2.8  n  N  {^y/y) 

Ey  =  rms  emittance  in  y-direction 
A/g  =  number  of  Betatron  oscillations. 


This  is  peak  gain  at  the  optimum  resonance  condition. 


A  gain  plot  can  be  generated  to  illustrate  the  effects  of  emittance 
and  energy  spread  for  a  given  undulator  and  electron  beam.  This  analysis,  valid 
only  in  weak  fields,  is  a  one-dimensional  theory  where  all  transverse  effects  have 
been  neglected.  However,  the  size  of  the  electron  beam  relative  to  the  optical 
mode  size  is  incorporated  by  including  the  filling  factor.  The  five  methods  used  are 
summarized  below  in  Table  II. 


TABLE  II. 


FEX 

WFS 

SIM 

INT 

DAT 

high  gain 

yes 

yes 

yes 

yes 

no 

mode  distortion 

yes 

yes 

no 

no 

no 

energy  spread 

yes 

no 

yes 

yes 

yes 

strong  fields 

yes 

yes 

yes 

no 

no 

harmonics 

yes 

yes 

yes 

yes 

yes 

The  high  gain  theories  are  non-linear  in  current  j  and  are 
determined  self-consistentty.  A  theory  with  mode  distortion  describes  a  self- 
consistent  calculation  of  the  transverse  optical  wavefront  without  the  use  of  the 
filling  factor.  Modeling  the  electron  beam  to  include  energy  spread  effects  is 


60 


incorporated  in  all  of  the  methods  except  WFS.  The  strong-field  theories,  FEX, 
WFS,  and  SIM  are  able  to  calculate  gain  in  strong  nonlinear  fields.  All  methods  of 
theoretical  analysis  include  harmonic  coupling  using  Bessel  functions  to  express 
the  reduced  coupling.  The  gain  for  each  theory  is  calculated  with  the  optimum 
resonance  condition. 

These  five  methods  each  have  their  own  merit  in  predicating  gain. 
The  FEX  code  incorporates  as  much  physics  as  is  feasible,  but  requires  a  CRAY 
computer  and  is  not  easy  to  use.  The  SIM,  INT  and  3-D  WFS  methods  require 
only  a  personal  computer,  but  do  not  include  the  full  range  effects  in  FEX.  On  the 
other  hand,  the  parameterized  gain  formula,  DAT,  is  valid  only  for  small-signal  gain 
in  weak  fields,  but  requires  only  a  calculator  for  evaluation. 

F.  DISCUSSION  OF  THE  PLOTS 

Plots  of  gain  for  the  three  micro-undulators  are  shown  in  Figures  1 5  through 
20.  In  Figure  1 5,  gain  is  plotted  against  the  rms  energy  spread  for  the  undulator 
using  the  fundamental.  The  five  methods  of  analysis,  FEX,  WFS,  SIM,  INT  and 
DAT  are  shown.  Only  a  single  WFS  data  point  is  shown,  indicated  by  an  X,  for  a 
monoenergetic  electron  beam.  The  emittance  is  e=>./4,  and  the  filling  factor  found 
to  be  9'  =0.25.  As  expected,  the  gain  decreases  for  all  four  methods,  FEX,  SIM, 
INT  and  DAT  when  the  energy  spread  is  increased  to  0.3%.  The  characteristic 
energy  spread,  shown  on  the  horizontal  axis  by  an  arrow,  indicates  where  gain 
should  be  seriously  degraded  by  the  energy  spread  effect  because  the  electron 
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GAIN  VS  ENERGY  SPREAD 


Figure  15.  Gain  Curve  for  Micro*undulator  Design  Operating  on  the 
Fundamentai  where  Energy  Spread  is  Varied. 

phase  is  roughly  equal  to  X.  FEX  and  WFS  have  nearly  equal  gain  values, 
because  they  represent  a  complete  self-consistent  FEL  oscillator  solution  where 
optical  distortion  of  the  optical  wavefront  is  taken  into  account  The  SIM  and  INT 
curves  are  below  FEX  for  the  entire  energy  spread  range,  because  they  do  not 
indude  mode  distortion.  The  DAT  theory  does  not  include  mode  distortion  or  self- 
consistent  high  gain,  and  therefore  lies  below  all  other  theories. 
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Gain  vs  optical  power  for  the  fundamental  undulator  with  nearly  perfect 
electron  beam  quality,  energy  spread  of  0%  and  e=A./4,  for  FEX,  SIM,  and  WFS 
are  shown  in  Figure  16.  The  gain  for  the  three  methods  starts  to  degrade  at 
approximately  the  same  optical  power.  The  characteristic  power  where  the 
beginning  of  saturation  should  occur  is  indicated  by  a  star.  The  FEX  and  WFS 
curves  are  above  the  SIM  curves,  because  mode  distortion  increases  coupling  at 
all  powers.  The  3-D  codes  agree  throughout  the  entire  power  range  except  when 
the  optical  field  becomes  very  strong.  There,  an  exact  calculation  of  gain  from  any 
FEL  model  becomes  difficult. 

Figure  1 7  displays  gain  vs  energy  spread  for  the  second  undulator  design 
which  operates  on  the  third  harmonic  with  design  parameters  given  in  Table  I.  The 
gain  curve  is  similar  to  the  previous  design  except  that  the  filling  factor  is  9'  =0.1 . 
Similar  trends  in  the  gain  degradation  and  saturation  are  seen  in  Figures  1 5  and 
1 6,  because  of  the  attributes  of  the  theoretical  models  used.  However,  this  design, 
operates  on  the  third  harimonic  with  only  A/=49  periods  compared  to  A/=413 
periods  for  the  fundamental  design.  The  sensitivity  of  gain  to  energy  spread  is 
reduced  significantly  as  shown  in  Figures  15  and  16.  The  gain  vs  optical  power 
plot,  shown  in  Figure  18,  illustrates  the  additional  advantage  of  using  harmonics 
since  the  characteristic  power  for  saturation  6x1 0®  W/cm^,  is  approximately  a  factor 
of  15  higher  than  in  the  fundamental  design.  As  seen  in  Figure  18  some 
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GAIN  VS  POWER 


INPUT  POWER  (W/cm  *  ) 


Figure  16.  Gain  Curve  for  Micro-undulator  Design  Operating  on  the 
Fundamental  where  Optical  Power  is  Varied. 
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OPTICAL  GAIN  (%) 


GAIN  VS  ENERGY  SPREAD 


1  I  ■  '  I  a'  I  "  •  ■!  *  —I  •  ■  I 

0.0  0.2  ^  0.4  0.0  0.0  1.0 

ENERGY  SPREAD  (%) 


FIguro  17.  Gain  Curva  for  MIcro-undulator  Design  Operating  on  the  Third 
Harmonic  whore  Energy  Spread  is  Varied. 
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deviations  between  FEX  and  WFS  occur  at  very  high  powers  where  calculations 
are  difficult. 

Extending  the  comparison  to  the  third  undulator  operating  in  the  ninth 
harmonic,  Figure  19  shows  the  same  general  features  as  in  the  previous  gain  vs 
energy  spread  plots.  Only  N=^9  periods  are  required  here  compared  to  N=A9 
periods  for  the  third  harmonic  design  and  A/=413  in  the  fundamental  design.  The 
gain  vs  optical  power  plot,  Figure  20,  illustrates  the  same  shape  for  the  three 
methods  as  the  third  harmonic  design.  However,  the  characteristic  power  for 
saturation  7x10®  W/cm^  is  approximately  a  factor  of  20  higher  than  found  in  the 
fundamental  design. 
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OPTICAL  GAIN  (%) 


GAIN  VS  ENERGY  SPREAD 


ENERGY  SPREAD  (%) 


Figure  19.  Gain  Curve  for  Micro-undulator  Design  Operating  on  the  Ninth 
Harmonic  where  Energy  Spread  is  Varied. 
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V.  CONCLUSIONS 


Preliminary  designs  of  Compact  Free  Electron  Lasers  are  analyzed.  This 
work  contributes  to  the  operational  advancement  of  Compact  FELs  by  illustrating 
the  advantages  of  using  harmonics  where  size  and  complexity  would  be 
substantially  reduced.  The  device  would  be  available  to  a  wider  range  of  users 
with  many  potential  applications  in  scientific  and  modern  technology.  Three  micro- 
undulator  designs  are  explored  to  reach  an  optical  wavelength  of  0.4  microns  with 
an  electron  beam  of  20  Mev.  The  use  of  harmonics  at  short  wavelengths  improves 
the  undulator  design  with  longer  periods  that  would  be  easier  to  fabricate.  The 
use  of  shorter  undulators  minimizes  the  penalty  from  emittance  and  energy  spread 
inherent  in  all  electron  beams.  For  fixed  optical  wavelength,  and  an  undulator 
using  has  monies,  the  electron  beam  energy  can  be  reduced  thereby  lowering  the 
radiation  shielding  requirements  for  personnel.  Operation  over  a  wider  range  of 
wavelengths  and  the  increase  in  optical  power  available  at  saturation  provide 
additional  benefits  of  using  harmonics.  The  value  of  using  harmonics  to  improve 
undulator  designs  are  evident  by  performance  calculations  conducted  by  various 
computer  codes,  and  should  be  pursed  to  extend  this  new  scientific  tool  for  future 
applications. 
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